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Background: Lutetium thermometry has been used to analyze Oklo natural nuclear reactor zones but leads to 
widely varying and puzzling predictions for the temperatures To which in turn impacts bounds on time variation 
of the fine structure constant a. 

Purpose: We revisit results for reactor zone RZIO in light of new measurements of the isomer branching ratio 
in ^^''Lu neutron capture at 5 and 25 keV. 

Method: We recalculate predictions for To as a function of using realistic models of the Oklo neutron flux. 

Results: We find To ~ 100 ±30 C using a new value of B^ , in contrast to 350 < To < 500 C using the evaluated 
value at thermal energy. 

Conclusions: Lutetium thermometry can be applicable to analyses of Oklo reactor data, but a better measure- 
ment of B^ with thermal neutrons is needed to confirm the reliability of temperature predictions. 
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I. PROBLEM OF TEMPERATURE IN THE OKLO REACTORS 



The Oklo natural reactors (see e.g. Naudet [T|) have proven to be one of the more sensitive terrestrial testing grounds 
for studying time variation of the fine structure constant a 0-01 or time variation of the dimensionless quark mass 
parameter Eg of QCD 8]. The position of the first neutron resonance in ^^^Sm strongly influences the effective 
cross section a for neutron capture leading to burnup of ^'^^Sm. The present day resonance energy Eo=97.3-meV is 
well known, but in the past may have been different, leading to a change (through a) in the isotopic abundances 
of various Sm isotopes in Oklo reactor wastes. From a modeling of the reactor parameters, and in particular from 
knowing the reactor temperature, thermal and epithermal neutron spectra can be derived which can in turn be used 
to predict whether the ^^^Sm cross section has changed since the reactors stopped operating. While the majority of 
Oklo analyses have been consistent with no change, a positive effect continues to be argued for from astronomical 
observations f£|. 

Oklo analyses all reveal that the bounds on variation of a depend significantly on the assumed reactor operating 
temperature Tq- In Refs. [l-lSl Maxwellian thermal neutron fluxes were used, while in Refs. 0-0] several realistic 
models of the reactors zones were elaborated with the help of modern neutron transport codes. For zone RZ2, with 
indirect arguments, Damour and Dyson Q allowed for a broad interval of Tq from 450 to 1000 C. Y. Fujii et al. [1] 
took the interval from 180 to 400 C for zones RZ3, RZIO, while Gould et al. @ preferred lower values from 200 to 
300 C for zone RZIO. Petrov et al. [6j argued for the value Tq = 452 ± 55 C as the temperature value at which their 
model of the active core RZ2 became critical Q. 

Three studies have attempted to bound temperatures more directly using ^^^Lu/^'''^Lu thermometry. Holliger and 
Devillers [l0|, gave temperatures Tq = 260 and 280 C for RZ2 and RZ3, and recently OneginH found Tq = 182 ±80 
C for RZ3. However, for RZIO, one of the most well characterized zones, Hidaka and Holliger[llj succeeded in getting 
a result for only one sample Tq — 380 C while for the other three samples only a (surprisingly high) lower bound. 
To > 1000 C, was obtained. Obviously, the temperature Tq of Oklo reactors remains a very uncertain parameter 



and the results of Ref. ll| are especially puzzling. In what follows we will discuss some possible ways to improve the 
situation. 

In Oklo studies, the effective capture cross section is introduced as 

n{v) a{v) V dv/nvQ. (1) 



Here n — n{v)dv is the total neutron density, and 1)0=2200 m/sec is the velocity of a neutron at thermal energy 
0.0253 eV. As an integrated quantity, a is not dependent on neutron velocity v (or neutron energy E) but if some 
nuclide, like ^''^Lu, has a neutron resonance close to thermal energy, ^{To) may depend on temperature Tq through a 
possible temperature dependence of the neutron density, n{E,To)- Examples of calculated curves n{E,To) for zones 
RZ2 and RZIO are given in Rcf. 5]. Also an effective neutron flux density is introduced as $ = n{v)vodv. This 
flux is different from the integrated 'true' flux $ = n{v)vdv, but the reaction rate R ^ = is the same 
since the average cross section a = n{v) (j{v) v dv / n{v) v dv. 

We will follow the description of the Lutetium thermometry given by Holliger and Devillers [l3l . The first suggestion 
to use Lutetium as a sensitive indicator of the temperature of neutron spectra in reactors was made apparently by 



Westcott [12 1 . 



II. i^'^LU/^^'^LU THERMOMETRY 

The rare earth element Lutetium has one stable isotope, ^^^Lu, with natural abundance 97.401 %, and and a second 
very long-lived isotope ^''^Lu with half-life ti/2 — 37.6 Gyr and present day natural abundance 2.599 % [l^. Neither 
isotope is produced in the fission of Uranium or Plutonium. 

^''^Lu has an exceptionally large thermal capture cross section, erg = 2090 ± 70 b [3l due to a strong resonance 
at 141 meV. This leads to a strong temperature dependence in the rate at which ^^^Lu burns up. Following neutron 



We note that the Tq parameter depends strongly on the H/U atomic ratio, which in their active core model had rather low value, 
H/U=5.9, as compared with H/U=15.6 in a subsequent model of Onegin Q. 
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capture, -'^''^Lu transforms into the stable nuclide -'^^^Hf after the /3 decay of the product -'^''^Lu . Neutron capture by 
^^^Lu has two branches with much smaller cross sections. The dominant branch with ct™ — 16.7 ± 0.4 b leads to 
an isomeric state ^^^™Lu which decays to ^^^Hf with half-life of 3.6 hr. The minor branch with erf = 6.6 ± 1.3 b [3| 
leads directly to the ground state of ^''^Lu. While the large cross section for neutron capture in ^^^Lu serves to deplete 
the small fraction of ^^^Lu, the minor capture branch on the much larger fraction of ^^^Lu works to restore ^^^Lu. 
Defining (75 = erf + cr™, the important parameter determining the balance between depletion and restitution of ^^^Lu 
is then — a^ja^^ the branching ratio parameter. According to these evaluated cross sections — 0.28 ± 0.05 
at thermal neutron energies. 

Introducing the atomic number densites Ni and the reaction rates — (the latter play roles analogous to the 
roles of decay constants in the radioactive decay) and neglecting for now the /3 decay of the 'stable' ^''^^Lu during 
the relatively much shorter time of Oklo reactor criticality, we write the coupled differential equations for the time 
evolution of the number densities Ni(t) of our two isotopes of interest: 
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dt 

dN5 



-XgNe + Bn^N^ (2) 
dt 

with initial conditions iV6(0) = TVg exp( '^'"^ ) and N^i^O) — N^. In these equations, the subscripts 5 and 6 refer to 

1/2 

^^^Lu and ^^^Lu, respectively, Nq and N^aie the present day natural abundances of lutetium isotopes and D = 2 
Gyr is the date of the Oklo phenomenon. 

This system of equations has an analytical solution which we write for the time duration ti of reactor criticality: 

iVel^i) =A^6°exp(^)exp(~A6ti) + iV0S9-A^[exp(-A5ti)-exp(-A6ti)] (4) 

tl/2 M — ^5 

7V5(ti) = 7V5"cxp(-A5ii). (5) 

The two terms in the right hand side of Eq.(4) represent the burnup of the initial ^''^Lu with cross section and its 
partial restitution after burnup of ^^^Lu with the partial cross section (t|. Taking the ratio and accounting for the /3 
decay of ^^^Lu after shut down of the reactor, we obtain the present day ratio of the lutetium isotopes in Oklo wastes 
as: 

now = -| exp - - (75 $ti + ^ \ 1 - exp(-0T6 - (75 $ti exp , 6 

which depends on the temperature Tq through the temperature dependence of (Tg. 

The neutron fluence $^1 is a well characterized parameter in Oklo studies. The effective cross sections 175 and erg 
are calculable using known resonance cross sections of Ref. 14] together with neutron densities within realistic reactor 
models. The calculation therefore leads to a clear prediction of the present isotopic ratio if the value of the branching 
ratio parameter B^ is known. 

The determination of B^ has been the subject of considerable experimental effort in connection with astrophysical 
studies of the s-process in stars. The difficulty in making a precise determination arises from the fact that its value is 
obtained by taking the difference between two cross sections, as and cr™, of similar magnitude. In a 1988 summary of 
measurements [la], a value anywhere in the range 0.16 to 0.38 seemed possible. To address the problem, two precision 
measurements of the "'^''^Lu total capture and partial activation cross sections at 5 and 25 keV were recently carried out 
[1IIT3. The conclusions of these studies were that [16] B^ = 0.117±0.046 at 5 keV and B^ = 0.143±0.020 at 25 keV. 
These values are about a factor of two less than the value for thermal neutrons, 0.28. Such a big difference looks rather 
puzzling even taking into consideration a partial contribution of p-wave neutrons, and raises the question of whether 
the thermal value is actually correct. There is some experimental evidence for a smaller value. A spectroscopic study 
[m of thermal neutron capture found B^ = 0.13 ± 0.03. Also, the only neutron transmission measurement for which 
there are data with an enriched ^^^Lu sample, (Baston et al. [l^), hints at a smaller value in the following way. 
The thermal capture cross section they determined, (75 = 23 ± 3 b, was based on subtracting a potential scattering 
contribution of 5.5 b from their measured total cross section. However, the potential scattering cross section is known 



^ Another parameter, the isomeric ratio, IR is also used in the literature: IR = 1 — . 

^ The age of the Oklo natural reactors is debated in the literature: Fujii et al 3) and Dyson and Damour use 2 Gyr. Naudet cites 
1.95 Gyr Q] and 1.8 Gyr is used in Ref. Q. 
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today to be much larger, (t„ = 7.2 ± 0.4 b [l^l- Using this new value leads to (75 = 21.3 b, and correspondingly a 
smaller Ba = 0.21 ± 0.11 b. 

While this latter result has very large uncertainty, we believe there is reason to think the value of for thermal 
neutrons is not that well known, and merits further study. In what follows we evaluate Oklo reactor temperature 
predictions for a range of values from 0.10 to 0.30, and then focus on the most precise value measured to date, the 25 
keV value cited above. 



III. To CALCULATION FOR ZONE RZIO 

To solve equation 6, we first evaluate the effective capture cross sections a for ^^^Lu and a for -'^''^Lu using the T = 
20 to 500 C neutron spectra derived from MCNP calculations in our earlier work Q . The procedure was identical to 
that outlined there except here the Lu-resonance shift due to a change in a is taken to be zero. For I'^^Lu we include 
all resonances up to 49.4 eV along with two sub-threshold resonances and the strong resonance at 96.69 eV. For ^^^Lu 
we include all resonances up to 52.13 eV. In practice, the contribution of the Eq = 143-meV resonance dominates for 
-^^^Lu, all other resonances contribute only a few percent. 

TABLE I. (Lutetium cross sections 0-5 and aa for the Oklo RZIO reactor at temperatures To from C to 600 C (see text). 



To (C) 


<T5 (kb) 


<T6 (kb) 





0.115 


4.216 


20 


0.115 


4.487 


100 


0.115 


5.359 


200 


0.115 


6.310 


300 


0.114 


7.013 


400 


0.114 


7.544 


500 


0.114 


7.715 


600 


0.114 


7.750 



TABLE II. RZIO borehole SF84 data on Lutetium isotopic abundances [Tl|] and neutron fluences from [^. For RZIO ti = 850 
kyr. The metasample data point is a weighted average of the data for the four individual samples. 



Sample 


Elemental Lu (ppm) 


N^lnow) 


Na (now) 


Nq {now) 
Nt^ (now) 


cj>ti (kb)-^ 


natural 




97.401 


2.599 


0.02668 




1469 


0.934 


99.695 


0.303 


0.00304 


0.475 


1480 


0.876 


99.479 


0.521 


0.00524 


0.915 


1485 


1.190 


99.866 


0.134 


0.00134 


0.645 


1492 


1.710 


99.377 


0.623 


0.00627 


0.585 


met a 








0.00418 


0.650 



The resulting cross sections are shown in Table I. Here we also include values for Tq — C and Tq — 600 C, derived 
from a power series extrapolation of the 20 C to 500 C results. 

The Lu isotopic abundance data for reactor zone RZIO are given by Hidaka and Holliger (HH) [11] and are shown 
in Table II for four different borehole depths. We see "'^''^Lu is strongly depleted in all samples although the ratio 
iV6/iV5(now) varies substantially, indicating, as noted by HH, heterogeneous operation of the reactor zone, or dif- 
ferences in isotope retention following reactor shut down. Lacking detailed further information on this situation, it 
appears more useful to work with a meta sample averaged over the four samples that are available for analysis. This 
approach was shown to be successful for Sm data [5^. Weighting each sample by its Lu elemental concentration, we 
accordingly find the meta sample ratio / N^{now) — 0.00418. 

This meta sample value can now be compared to the result of evaluating equation 6 for a range of B^ values and 
reactor temperatures. A plot of this comparison is shown in Fig. 1 for B^ values 0.1, 0.15, 0.2, 0.25 and 0.3, and 
temperatures to 600 C. Of interest is that for a given B^ only a small range of isotope ratio values actually lead 
to any solution at all. Below about 350 C, the temperature is well constrained by the isotopic ratio data. At higher 
temperatures, however, the curves flatten out, leading to much reduced sensitivity, and above 500 C there is actually 
no solution. 

The intersection of the meta sample isotope ratio value with each curve yields a temperature prediction for the 
RZIO reactor. A plot of these intersection values is given in Fig. 2 from which a temperature Tq for RZIO can be 
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FIG. 1. (Color online.) Plot of solutions of Eq. 6 as a function of temperature To for values 0.1, 0.15, 0.2, 0.25 and 0.3, 
and temperatures to 600 C. The meta sample value is the horizontal line Ne/N3{now) — 0.00419. 
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FIG. 2. (Color online.) Plot of the intersection temperatures for the meta sample as a function of the isomer branching ratio 
. Note there is no solution above 500 C. 



read off if is known. 
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As noted earlier, the evaluated thermal neutron energy value [11] is = 0.28±0.05 which (barely) yields a solution 
350 <To < 500 C. The most accurately determined (£;„ = 25 kcV) value of 0.157 ± 0.023 from the study [13] 
gives a quite different value Tq — 100 ± 30 C. 



IV. CONCLUSIONS 



In this work we have identified the branching ratio parameter B^ as a possible source of problems in extracting 
reliable temperature bounds from Lutetium thermometry for the Oklo natural nuclear reactors. Our preference is for 
a lower value, but clearly there is impetus for additional study to constrain Tq better. 

A first possibility is to carry out an improved measurement of B^ at thermal energies. A ten percent measurement 
of B^ is desirable and will certainly require carrying out a high precision total cross section measurement. An accuracy 
of 1% has been achieved at RPI [2l[ for other rare earth element total cross sections, and similar accuracy should be 
feasible with an enriched ^^^Lu sample. The ^^^Lu isomer activation cross section is less critical but should also be 
improved bearing in mind that the differences in the existing experimental values considerably exceed the reported 
uncertainties. 

There remains of course the issue of whether the meta sample, even though averaged over the samples of the RZIO 
active core, still deviates from the true value 2 Gyr ago. Post Oklo migrations of elements in and out of the core have 
been discussed in the literature [2^ . The apparent strong retention of other rare earth samples in the active cores 
would seem to argue against selective migration of Lutetium alone, but additional samples and data would clearly be 
valuable in clarifying the migration issue. 

Finally there is the possibility that post processing of the Lutetium isotopes has occurred due to other nuclear 
reaction channels. This has been a topic of on-going interest in astrophysics, where an excess of ^^^Hf in older 



meteorites had been a long standing puzzle. See Thrane et al. |23| for a summary of the current situation. It is now 
generally agreed that the excess is not due to an incorrect value for the ^ ''^Lu lifetime, but is instead due to excitation 
of the ^^^""Lu isomer by gamma radiation. The source of the gamma ray flux is debated, but as discussed in Ref. 
[2^, a cosmic ray spray from a nearby supernova about 5 Gyr ago is a plausible mechanism. The fluorescence cross 
section for ^^^Lu isomer production has been intensively studied by a number of groups (see Mohr et al. [13]). Cross 
section results in the literature vary by many orders of magnitude depending on the energy of the gamma radiation. 
Although the gamma-ray fluxes in nuclear reactors are extremely low as compared to those in supernovae, additional 
measurements of the ^^^Lu photoactivation cross sections and calculations of gamma-ray fluxes in Oklo reactors are 
still desirable for a better understanding of the Oklo phenomenon. 
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